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tice, using a renormalization group improved gauge action and the clover improved Wilson quark 
action. In the simulations along lines of constant mps/my, we calculate the Taylor expansion 
coefficients of the heavy-quark free energy with respect to the quark chemical potential (/i^) up to 
the second order. By comparing the expansion coefficients of the free energies between quark(2) 
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1. Introduction 

We study QCD thermodynamics at finite temperature (r) and quark ciiemical potential (Hq) 
with two-flavors of dynamical quarks. Then, we calculate the Taylor expansion coefficients of 
physical quantities in the simulations at = 0, and investigate thermodynamic properties at small 
Hq region. In this report, we present current status of two topics: heavy-quark free energy and 
fluctuation at finite Hq. The former is related to the inter-quark interaction in quark-gluon plasma, 
and the latter is related to the existence of the critical point in {T^jXq) plane. 

2. Lattice action and simulation parameters 

We employ a renormalization group improved gauge action and a clover improved Wilson 
quark action with two flavors. The numerical simulation is performed on a lattice with a size of 
X Nt = Ky' along lines of constant physics, i.e. lines of constant mps/my (the ratio of 
pseudoscalar and vector meson masses) at T = in the space of simulation parameters. Two values 
of mps/my are taken: 0.65 and 0.80 with the temperature range of T /Tpc ~ 0.82-4.0 and 0.76-3.0, 
respectively, where Tpc is the pseudo-critical temperature along the line of constant physics. The 
number of trajectories for each run after thermalization is 5000-6000, and we measure physical 
quantities at every 10 trajectories. Details of the lines of constant physics with the same actions are 
summarized in Ref. [jl], ^. To calculate derivatives of the quark determinant with respect to pLq, we 
use the random noise method introduced in Ref. ^ with the number of noise of 100-200. 

3. Heavy quark free energy and Debye screening mass 

The heavy-quark free energy in the QCD medium is one of the important quantities to char- 
acterize the properties of the quark-gluon plasma. Especially, the properties of the free energies 
at finite T and jiq may be related to the fate of the charmoniums and bottomoniums in relativistic 
heavy ion collisions. Precise studies of the heavy-quark free energy in two-flavors QCD have been 
done with the improved Wilson quark action at /i^ = [Q]. Properties of the free energy between 
static-quark {Q) and -antiquark {Q) at finite /i^ have been previously investigated by using the Tay- 
lor expansion method using an improved staggered quark action [^. In this section, we show the 
Taylor expansion coefficients of the free energy between not only Q and Q, but also Q and Q, up to 
2nd-order of }JLq. 

The free energy of static quarks on the lattice is described by the correlations of the Polyakov 
loop: il(x) = n^Li Uis,{x,ii) where the ?7^(t,x) G SU(3) is the link variable. With an appropriate 
gauge fixing (e.g. the Coulomb gauge fixing), one can define the free energy in various color 
channels separately: the color singlet QQ channel (1), the color octet QQ channel (8), the color 
anti-triplet QQ channel (3*), and the color sextet QQ channel (6), given as follows. 

^-F\r.T)lT 



= i(Tra^(x)a(y)), (3.1) 
= i(Trn^(x)Trn(y))-l(Tra1^(x)a(y)), (3.2) 
= -i-(Tra(x)Tra(y)) + -i-(Tra(x)a(y)), (3.3) 
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Figure 1: Results of Vq (left) and V2 (right) for QQ channel above Tpc at mps /my = 0.80. 



^-F^\r,T)/T ^ LTrQ(x)Trn(y)) - l(Tril(x)Q(y)), (3.4) 

where r = |x — y|. 

For T > Tpc, we introduce normalized free energies {V^ ) such that they vanish at 
large distances. This is equivalent to define the free energies by dividing the right-hand side of 



Eq. (3J)-(3^) by (Tr^l)(Tr^l^) for QQ free energies and (TiQ.)^ for QQ free energies. 

The Taylor expansion of normalized free energies with respect to Hq/T is described as a power 
series up to 2nd order: 

r,M,) =v^ + vf (^) + vf (^)' + (3.5) 

where M is the color channel. The coefficients, v^, can be evaluated by expanding the quark 
determinant of partition function in powers of /i^, then the normalization of the free energies by 
(TrH) is also taken order by order of /i^. 

One should note that the color singlet and octet channels do not have the odd orders in the 
Taylor expansion since the free energies for both channels are symmetric with respect to jj.^. In 
other words, the free energies between Q and Q are invariant under the charge conjugation. On the 
other hand, the color sextet and antitriplet channels has the odd orders since the QQ free energies 
are not invariant under the charge conjugation. 

The expansion coefficients of the normalized free energies for the color singlet and octet QQ 
channels are shown in Fig. [l| for (left) and vf (right) at mps/my = 0.80 for several temperatures. 
Those for the color sextet and antitriplet QQ channels are shown in Fig. § for (left) and (right) 
and in Fig. || for vf. It have been found in Ref. that the inter-quark interaction is "attractive" 
in the color singlet and antitriplet channels and is "repulsive" in the color octet and sextet channels 
at = 0. We find that, both in high and low temperatures, the sign of vf is the same with that 
of v^, whereas the sign of a is the opposite of that of Vq, i.e. v^^ • > (only for QQ free 
energies) and vf - v^ <0. This means that the inter-quark interaction between Q and Q becomes 
weak, whereas that between Q and Q becomes strong in the leading-order of Hq. In other words, 
QQ (QQ) free energies are screened (anti-screened) by contributions of the internal quarks induced 
by finite Hq. 
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Figure 2: Results of v^ (left) and vf (right) for QQ channel above Tpc at mps /my = 0.80. 
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Figure 3: Results of v^ for QQ channel above Tpc at mps /my = 0.80. 

In order to study the screening effect in each color channel, we fit the normalized free energies 
by a screened Coulomb form, 



V''{r,T,pi,)=C{M) 



-mD(T,llq)r 



(3.6) 



where aeff(7",/x^) and mQ{T,ji.q) are the effective running coupling and the Debye screening mass, 
respectively. We assume that contributions of finite /i^ appear only in ttgff and nio. The Casimir 
factor C{M) = (L^^i t" ■ ff)^ for color channel M is given by 

4 , - 1 . . 1 , 2 

~- C(6) = -, C(3*) = --. (3.7) 



CI 



C8 



3' ' ' 6' ' ' 3' ' ' 3 
We assume that mB{T,}Xq) is also expressed as a power series of IXq/T: 



(3.8) 



where we use the fact that the Debye screening mass does not have the odd powers in the Taylor 
expansion because it is related to the self-energy of the two-point correlation of the gauge fields 
which is symmetric when Hq —piq- Relations of coefficients between V'^ and mr, are given by 
comparing each order of jlq. At jiq = 0, the relation is the same as that adopted in Ref. [^: 



vo(r,r) = C(M)^^5^^'-'"^»(^)'-. 



(3.9) 
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Figure 4: Results of mo o (left) and mo 2 for each color channel at mps/my — 0.65. Dashed lines are 
prediction from a leading-order of the thermal perturbation theory with the renormalization points of K — nT, 
2kT and3;rr. 



For the second order coefficients, we obtain, 



V2 

— ~ -niDir. 



(3.10) 



At large distances, we estimate the coefficients of Debye mass by fitting the normalized free ener 
gies for each color channel with the formulae (I 



and ( P.IOD . 

Figure ^ shows the results of the mj)fl{T) (left) and mD,2{T) (right) for mps/my = 0.65. Sim- 
ilar behavior in both results are obtained for nips/niy = 0.80. We find that there is no significant 
channel dependence in both coefficients at sufficiently high temperatures {T > 2Tpc). In other 
words, the channel dependence of the free energy at high temperature can be well absorbed in the 
Casimir factor. 

Let us compare the Debye screening mass at finite /i^ on the lattice with that predicted in the 
thermal perturbation theory. The 2-loop running coupling is given by 



jSoln 



mv^ininm 



A/ 



(3.11) 



where K is the renormalization point. The argument in the logarithms can be written as k/A = 
{k/T)[T /Tp,){TpjA) with A = A^^"^ ~ 261 MeV ||] and Tpc ~ 171 MeV [§. We assume k to 
be in a range TiT to SttJ. Therefore, g2\ can be viewed as a function of T /Tp^-. In the leading order 
of the thermal perturbation theory, the Debye screening mass with g2\ is given by 



1 /2 

m)P{T,H,)=g2m[{l + ^)T' + ^,nl] ■ 



(3.12) 



Therefore, the leading-order expansion coefficients in the thermal perturbation theory are given by 



m 



LO , 
Dfi\ 



{T) = \l + ^g2i{T)T 



Nf 



4712 



82i{T)T. 



(3.13) 



The dashed lines in Fig. ^ are the results of and for K = nT, 2nT and 3nT. We find 
that the lattice results are larger than the leading-order thermal perturbation for both coefficients 
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of the Debye mass. Since it is known that the next-to-leading-order in the thermal perturbation 
with respect to T compensates for such discrepancy in the case at /x^, = the higher order 
contributions in the perturbation theory with respect to Hq could help us to understand our results 
on the lattice. 

4. Hadronic fluctuations at finite fiq 

Hadronic fluctuations at finite density are observables closely related to the critical point in 
the [T^pLq) plane and may be experimentally detected by an event-by-event analysis of heavy ion 
collisions. The fluctuations can also be studied by numerical simulations of lattice QCD calculating 
the quark number and isospin susceptibilities, Xq and Xj. They correspond to the second derivatives 
of the pressure with respect to jiq and jii, where pLj is the isospin chemical potential. From a 
phenomenological argument in the sigma model, Xq is singular at the critical point, whereas Xi 
shows no singularity there 

In this section, we calculate the Xq and Xi and their second derivatives with respect to jXq and jU/ 
at iJ.q = IJ.I = 0. (Note that the odd derivatives are zero at {Xq = 0.) The details of the calculations are 
reported in Ref . [||] . The left panel of Fig. |5| shows Xqf^^ (circle) and Xi/T^ (square) at mps / my = 
0.8 and /i^ = /i/ = as functions of T /Tp^. We find that Xqf^^ and Xif^^ increase sharply at Tp^, in 
accordance with the expectation that the fluctuations in the QGP phase are much larger than those 
in the hadron phase. Their second derivatives d'^{Xq/T^)/d{iJ.q/T)^ and d^{xi/T'^)/d{iJ.q/T)'^ are 
shown in Fig. ^ (right). We find that basic features are quite similar to those found previously 
with the p4-improved staggered fermions [^. d'^{xi /T^) / d{pLq/TY remains small around Tpc, 
suggesting that there are no singularities in Xi at non-zero density. On the other hand, we expect 
a large enhancement in the quark number fluctuations near Tpc as approaching the critical point 
in the iT^jXq) plane. The dashed line in Fig. ^ (right) is a prediction from the hadron resonance 
gas model, d^Xq/^l^q ~ ^Xq/T^- Although current statistical errors in Fig. || (right) are still large, 
we find that d'^{Xq/T'^) / d{jjLq/TY near Tp^ is much larger than that at high temperature. At the 
right end of the figure, values of free quark-gluon gas (Stefan-Boltzmann gas) for Nt = 4 and for 
Nt = °° limit are shown. Since the lattice discretization error in the equation of state is known to be 
large at Nt =4 with our quark action, we need to extend our study to larger A^f for the continuum 
extrapolation. 

5. Summary 

We presented current status of thermodynamics of two-flavor QCD with the renormalization 
group improved gauge action and the clover improved Wilson quark action. Simulations were 
performed on a 16^ x 4 lattice and along lines of constant mps/my = 0.65 and 0.80. 

The properties of the heavy-quark free energies at finite jXq were studied in the Taylor expan- 
sion method up to 2nd order of Hq. We find that there is a characteristic difference between QQ and 
QQ free energies due to the first order coefficient of the Taylor expansion. It suggests that the inter- 
quark interaction between Q and Q (Q and Q) become week (strong) in the leading-order of Hq. We 
also extract the expansion coefficients of the Debye screening mass for each color channel up to 
2nd order of Hq. The Debye mass shows no significant channel dependence at T ^ 2Tpc, whereas. 



6 



Thermodynamics and heavy-quark free energies with improved Wilson quarks 



Y. Maezawa 



8 - , 

• 3Cq/7^ mpj,/my=0.80 

6 - ' 


SI! (N =4)- 


20 




ll ' 1 ' 1 ' 1 

■ 1 ■ d^iXJI^Wli /Tf 

■ vAij' /' ^"q ' 




■ • 






1 mp<,/my=0.80 


■ 

2 - • 


SB limit 


10 






^ ^ SB(N=4) 

■ - • * • 

■ ■ ■ ■ ■ ■ ■ ' 


■ 








fl 


SB limit 





















. 1 1 1 < 1 



1.5 

T/T 



2.5 0.5 



pc 



1.5 



Figure 5: Left: Quark number (circle) and isospin (square) susceptibilities at fi^ = fii = 0. Right: The 
second derivatives of these susceptibilities. 

we find disagreement witli leading-order predictions of the thermal perturbation theory. Since it is 
known that the next-to-leading-order of the thermal perturbation with respect to T well reproduces 
the Debye mass on a lattice at jx^ = [Q], the calculations of the higher order perturbation with 
respect to /i^ will help us to understand the results on the lattice. 

The fluctuations of quark number and isospin densities were also discussed. Although the 
statistical errors are still large, we find that Xq seems to increase rapidly near Tpc as Hq increases, 
whereas the increase of Xi is not large near Tpc. These behaviors qualitatively agree with the 
previous results obtained with the p4-improved staggered fermions. 
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